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A.  Introduction 


The  Loran  system  is  a  radio  aid  to  navigation  which 
utilizes  the  principle  of  hyperbolic  fixing.  The  locus  of 
points  for  which  the  difference  in  arrival  time  of  synchron¬ 
ized  signals  from  a  pair  of  transmitters  is  constant  det¬ 
ermines  a  hyperbolic  line  of  positions  (LOP) .  The 
intersection  of  two  hyperbolic  lines  of  position  from  two 
pairs  of  transmitters  (stations)  determines  a  position 
estimate  (a  fix) .  That  two  lines  of  position  and  hence 
two  pairs  of  stations  are  required  for  a  fix  does  not 
necessarily  mean  that  four  separated  stations  are  required, 
since  one  station  of  one  pair  may  be  colocated  with  one 
station  of  the  other  pair.  In  this  case,  a  Loran  triplet 
is  formed  (Figure  1) .  Triplets  may  be  joined  "end-to-end" 
by  station  colocation  to  form  a  Loran  chain  (Figure  2) . 

Loran  chains  are  common  on  both  the  East  and  West  Coasts 
of  the  North  American  continent. 

"Standard  Loran"  or  "Loran-A"  operates  at  a  fre¬ 
quency  just  below  2  MHz  and  is  now  only  used  in  the 
Pacific  area.  "Loran-C"  operates  at  100-kHz  and  is 
used  both  by  ships  and  aircraft  world  wide.  The  compu¬ 
tational  algorithm  and  programs  described  herein  can  be 
used  for  position  determination  with  Loran-C  triplets. 
Further  information  on  the  history,  development  and 
operation  of  the  Loran  systems  may  be  found  in  References 
1  and  2 . 
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(a)  Colocated  Master  Stations 


S 


1' 


M, 


(b)  Colocated  Slave  Stations 

M1 


(c)  Colocated  Master  and  Slave 


Figure  1.  Loran  Triplets 


B.  Program  Description 


The  Position  Determination  Program  requires  opera¬ 
tional  data  cards  (described  below)  for  on-location  position 
determination  from  Loran  triplet  time-difference  measure¬ 
ments.  Two  programs  are  used  to  prepare  an  operational 
data  card.  Operational  data  cards  should  be  prepared  and 
validated  prior  to  on-location  navigational  use.  Then, 
although  three  programs  are  used  in  total,  only  one 
program  is  required  for  position  determination,  the 
remaining  two  programs  are  used  to  prepare  operational 
data  cards . 

The  programs  are : 

Master  Data  Card  Program:  This  program  is  used  to  prepare 
master  data  cards.  A  master  data  card  requires  the  follow¬ 
ing  information  for  a  master  (M)  station/slave  (S)  station 
pair  (a  list  of  current  station  parameters  is  in  the 
appendix) s 

1.  A  M/S  pair  identification  number. 

2.  The  quantity  At  which  is  the  sum  of  the  coding 

delay  plus  the  one  way  base  line  time  in  microseconds. 

3.  The  latitude  and  longitude  of  the  master  station. 

4.  The  latitude  and  longitude  of  the  slave  station. 

Some  preprocessing  of  these  data  is  performed  before  the  master 
data  card  is  generated.  The  data  generated  require  only  one 
side  of  a  TI-59  magnetic  card  for  each  M/S  pair,  thus  a 
second  M/S  pair  may  be  placed  on  side  2  of  a  master  data 

card.  A  master  data  card  should  be  prepared  in  advance  for 
each  M/S  pair  that  might  be  used  within  an  area  of  operation. 


Operational  Data  Card  Program:  This  program  is  used  to 
prepare  an  operational  data  card  for  each  Loran  triplet 
that  might  be  used  within  an  area  of  operation.  Each 
operational  data  card  contains  data  merged  from  two 
master  data  cards  that  contain  the  M/S  pair  information 
for  each  pair  of  a  triplet.  These  merged  data  are  validity 
checked,  colocation  of  master  or  slave  determined  and 
encoded . 

The  inputs  required  to  generate  an  operational  data 
card  for  a  Loran  triplet  using  this  program  are  those 
contained  in  the  two  master  data  cards  that  are  associated 
with  the  triplet.  It  is  possible  to  prepare  and  store 
operational  data  cards  rather  than  master  data  cards.  This 
may  be  desirable  if  there  is  no  scarcity  of  cards  and  storage 
space,  however  the  number  of  possible  Loran  triplets  is 
considerably  larger  than  the  number  of  M/S  pairs. 

Position  Determination  Program:  This  program  is  used  for 
position  determination.  The  required  inputs  for  the  program 
are  the  indicated  time  differences  T  for  the  two  M/S 
pairs  of  the  triplet  and  the  data  contained  in  the  opera¬ 
tional  data  card  associated  with  the  triplet.  The  outputs 
of  the  program  are  the  computed  latitude  and  longitude  of 
the  Loran  fix.  Note:  Every  Loran  determination  gives 
two  solutions:  the  proper  solution  and  the  spurious  solution. 
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The  spurious  solution  can  almost  always  be  rejected  by 
inspection.  However,  if  the  stations  of  the  Loran  triplet 
are  nearly  aligned  this  may  not  be  the  case,  and  additional 
information  will  be  required  to  determine  the  proper 

solution. 


I- 


C.  TI-59  Calculator  Programs 
1.  User  Instructions 


Master  Data  Card  Program  (MDCP) 


Input  Output 


Step 

Instructions 

Data/Units 

Keys 

Data/Units 

0. 

Turn  the  calculator  on.  If  the  cal¬ 
culator  has  been  on,  turn  it  off, 
then  on  to  initialize  it. 

1. 

Read  side  1  of  the  MDCP  program  card. 

2. 

Press  CLR  and  then  read  side  2  of  the 
MDCP  program. 

CLR 

3. 

Input  a  unique  ID  number  for  the  loran 
pair.  See  Note  a. 

ID 

2nd, A' 

ID 

4. 

Input  the  coding  delay  At. 

At 

2nd,B' 

At 

5. 

Input  the  master  station  latitude. 

See  Note  b. 

*M 

A 

0 

6. 

Input  the  master  station  longitude. 

See  Note  b. 

XM 

B 

0 

7. 

Input  the  slave  station  latitude. 

*S 

C 

0 

8. 

Input  the  slave  station  longitude. 

Xs 

D 

0 

9. 

Run.  When  the  calculations  are 

E 

complete,  4  will  appear  in  the 
display;  then  the  display  will  blank. 

After  the  display  has  blanked,  enter 
a  card  in  the  read/write  slot  to 
record  the  master  data  card  for  the 
Loran  pair.  When  the  recording  is 
complete,  4  will  reappear  in  the 
display.  See  Note  c. 

Notes:  a.  Loran  pairs  are  coded  on  the  navigation  maps  using 
designators  such  as  9930X,  9930Y,  9930Z  and  9930W. 

It  is  suggested  that  the  ID's  for  these  pairs  be 
coded  as  9930.1,  9930.2,  9930.3  and  9930.4  respectively. 
However,  any  consistent  scheme  is  acceptable. 
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b.  The  format  for  latitude  and  longitude  data  input  is 
of  the  form:  DDD.MMSSFF,  where 

DDD  denotes  degrees 
MMM  denotes  minutes 
SS  denotes  seconds 
FF  denotes  hundredths  of  a  second. 

A  minus  sign  (-)  denotes  South  latitude  and  West 
longitude.  See  the  Appendix  with  regard  to  the 
latter  sign  convention. 

c.  When  the  display  has  blanked  after  4  has  been  displayed, 
the  calculator  is  prepared  to  write  a  card  side  with 
the  information  (data)  contained  in  bank  4.  When  a 
card  is  inserted,  it  will  be  recorded.  Do  not  restrict 
the  cards  motion  once  the  card  has  been  grasped  by  the 
read/write  mechanism.  When  recording  is  complete,  the 
card  will  stop  and  4  will  reappear  in  the  display. 

The  card  can  then  be  removed. 


Operational  Data  Card  Program  (ODCP) 


Step 


Instructions 


Input  Output 

Data/Units  Keys  Data/Units 


0.  Turn  the  calculator  on.  If  the  calculator 
has  been  on,  turn  it  off,  then  on  to 
initialize  it. 

1.  Read  side  1  of  the  ODCP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the  ODCP 
program  card. 

3.  Press  CLR  and  then  read  a  master  data  card 
containing  the  first  pair  of  a  Ioran  triplet. 

4.  Press  A. 

5.  Press  CLR  and  then  read  a  master  data  card 
containing  the  second  pair  of  a  Ioran 
triplet. 

6.  Press  B.  See  Note  a. 

7.  Run.  When  the  calculations  are  complete,  4 
will  appear  in  the  display;  then  the  display 
will  blank.  After  the  display  has  blanked, 
enter  a  card  in  the  read/write  slot  to 
record  an  operational  data  card.  When  re¬ 
cording  is  complete,  4  will  reappear  in  the 
display.  See  Note  b. 

Notes:  a.  Should  the  display  flash,  then  the  two  master  data 
cards  do  not  correspond  to  a  Loran  triplet.  Both 
the  latitude  and  longitude  of  the  colocated  stations 
must  be  identical  on  both  master  data  cards  in  order 
to  successfully  produce  an  operational  data  card. 

b.  When  the  display  has  blanked  after  4  has  been  displayed, 
the  calculator  is  prepared  to  write  a  card  side  with 
the  information  (data)  contained  in  bank  4;  and,  when 
a  card  is  inserted,  it  will  be  recorded.  Do  not 
restrict  the  cards  motion  once  the  card  has  been 
grasped  by  the  read/write  mechanism.  When  recording 
is  complete,  the  card  will  stop  and  4  will  reappear 
in  the  display.  The  card  can  then  be  removed. 


CLR 

CLR 

A 

cm 

B 

E 
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Position  Determination  Program  (PDP) 


Step 


Instructions 


Input  Output 

Data/Units  Keys  Data/Units 


0. 


1. 


2. 


3. 


4. 


5. 


6. 

7. 


8. 


9. 


10. 


Turn  the  calculator  on.  If  the  cal¬ 
culator  has  been  on,  then  turn  it  off, 
then  on  to  initialize  it. 

Read  side  1  of  the  PDP  program  card. 

Press  CLR  and  then  read  side  2  of  the  CLR 

PDP  program  card. 

Press  C  to  initialize,  and  then  read  the  C 

operational  data  card  for  a  Loran  tri¬ 
plet  that  you  are  receiving.  See  Note 
a  below. 

Input  the  observed  time  delay  from  the  T  A  T. 

first  Loran  pair.  A 

Input  the  observed  time  delay  from  the  B  T 

second  Loran  pair. 

Compute  the  1st  solution  latitude. 

Compute  the  1st  solution  longitude. 

Compute  the  2nd  solution  latitude. 

Compute  the  2nd  solution  longitude. 

See  Note  b. 

To  compute  new  solutions  for  new  time 
delays,  go  to  step  4.  To  compute  new 
solutions  for  a  new  triplet,  go  to 
step  3. 

Note  a.  If  the  data  card  does  not  read  (a  flashing  display  results) , 
press  CLR  and  repeat  the  step.  If  the  card  does  not  read 
on  the  2nd  try,  the  card  should  be  cleaned  (see  the  TI-59 
owners  manual)  or  a  new  data  card  prepared. 

Note  b.  Once  the  proper  solution  (solution  1  or  solution  2)  has 
been  determined  for  a  triplet,  then  only  step  6  (for 
solution  1)  or  step  8  (for  solution  2)  need  be  repeated 
after  inputing  new  time  delays  in  step  4  and  5. 


D  latitude 

R/S  Longitude 

E  Latitude 

R/S  Longitude 
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2.  Example  Problem  for  Loran-C  Pairs  9930X  and  9930Y 


An  example  use  of  the  Master  Data  Card  Program  using  coding  delays 
and  latitude  and  longitude  values  from  the  Appendix: 

Input  Output 

Step  Instructions  Data/Units  Keys  Data/Units 


0.  Turn  the  calculator  on.  If  the  cal¬ 

culator  has  been  on,  turn  it  off,  then 
on.  This  assures  that  it  will  be  in 
the  correct  initial  state . 

1.  Read  side  1  of  <he  MDC  program  card. 

2.  Press  CLR,  then  read  side  2  of  the  MDC 


program  card. 

CIA 

3. 

Input  the  ID  for  9930X.  Next  press 
2nd,  then  A. 

9930.1 

2nd, A' 

9930.1 

4. 

Input  the  coding  delay  At  for 

9930X.  Next  press  2nd,  then  B. 

36389.66 

2nd,B' 

36389.66 

5. 

Input  the  master  station  latitude, 
then  press  A. 

34.034604 

A 

0 

6. 

Input  the  master  station  longitude 
and  press  +/-  for  West,  then  press 

B. 

-77.544676 

B 

0 

7. 

Input  the  slave  station  latitude, 
then  press  C. 

46.463218 

C 

0 

8. 

Input  the  slave  station  longitude 
and  press  +/-  for  West,  then  press 

D. 

-53.102816 

D 

0 

9. 

Press  E.  After  4  has  appeared  in 
the  display,  enter  an  available 

E 

card  side  in  the  read/write  slot  of 
the  calculator.  After  recording  is 
complete,  4  will  reappear.  Then  re¬ 
move  the  card  and  label  the  recorded 
side  "9300X  MASTER." 

10.  Input  the  ID  for  9930Y,  next  press  9930.2  2nd, A'  9930.2 

2nd  and  then  A. 


Step 


Instructions 


Input  Output 

Data/Units  Keys  Data/Units 


11.  Input  the  codinq  delay  At  for 
9330Y,  next  press  2nd  and  then  B. 

12.  Input  the  master  station  latitude, 
then  press  A. 

13.  Input  the  master  station  longitude 
and  press  +/-  for  West,  then  press 

B. 


52541.31  2nd,B 1  52541.31 

34.034604  A  0 


-77.544676  B  0 


14. 

Input  the  slave  station 
then  press  C. 

latitude , 

41.151193 

C 

15. 

Input  the  slave  station 
and  press  +/-  for  West, 
D. 

longitude , 
then  press 

-69.583909 

D 

0 


0 


16.  Press  E.  After  4  has  appeared  in  E  0 

the  display,  enter  an  available 
card  side  (for  example  the  second 
side  of  the  card  used  in  Step  9)  in 
the  read/write  slot  of  the  cal¬ 
culator.  After  recording  is 
complete,  4  will  reappear  in  the 
display.  Then  remove  the  card  and 
label  the  recorded  side  "9930Y 
MASTER."  The  recorded  data  cards  (or 
card)  will  be  used  in  the  next 
example . 
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An  example  use  of  the  Operational  Data  Card  Program 


Instructions 


Input 

Data/Units 


Turn  the  calculator  on.  if  the  cal¬ 
culator  has  been  on,  turn  it  off,  then 
on  to  initialize  it. 

Read  side  1  of  Hie  ODCP  program  card. 

Press  CLR  and  then  read  side  2  of  the 
ODCP  Program  Card. 

Press  CLR  and  then  read  the  master 
data  card  for  station  9330X. 

Press  A. 

Press  CLR  and  then  read  the  master 
data  card  for  station  9930Y. 


Press  B.  See  Note  a  of  the  ODCP 
instructions  on  page  9. 

Press  E.  When  the  display  has 
blanked  after  4  has  appeared,  insert 
an  available  card  side  in  the  read/ 
write  slot.  After  recording  is 
complete,  4  will  reappear  in  the 
display.  Label  the  recorded  card 
side  "9330V9930Y  OPERATIONAL  DATA. 
Then,  for  identification,  label  the 
A  key  position  on  the  card  "9330X" 
and  the  B  key  position  "9330Y". 

This  card  will  be  used  as  the  data 
card  in  the  Position  Determination 
Program  example.  The  key  labeling 
relates  to  that  program. 


An  example  use  of  the  Position  Determination  Program 


Step 


Instructions 


Input 

Data/Units  Keys 


You  are  in  the  North  Atlantic  and  you 
are  receiving  9930X  and  9930Y  and  you 
wish  to  obtain  a  fix. 

0.  Turn  the  calculator  on.  If  the  cal¬ 
culator  has  been  on,  then  turn  it 
off,  then  on  to  initialize  it. 

1.  Read  side  1  of  the  PDP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the  CLR 

PDP  program  card. 

3.  Press  C  to  initialize,  and  then  read  C 

the  operational  data  for  the  triplet 

9930X/9930Y  from  the  card  side  re¬ 
corded  in  the  ODCP  example. 

4.  The  indicated  time  delay  is  28800  ys 
from  9930X. 

Input  the  indicated  time  delay,  then  28800  A 

press  A. 

5.  The  indicated  time  delay  is  49400  ys 
from  9930Y. 

Input  the  indicated  time  delay,  then  49400  B 

press  B. 

6.  Press  D  to  obtain  the  1st  solution  D 

latitude . 

7.  Press  R/S  to  obtain  the  1st  solution  R/S 

longitude. 

The  first  solution  is:  42°44'57"N  Latitude 

41°07'32"W  Longitude 

8.  Press  E  to  obtain  the  2nd  solution  E 

latitude . 

9.  Press  R/S  to  obtain  the  2nd  solution  R/S 

longitude . 

The  2nd  solution  is:  27°00'07"S  Latitude 

102 °27 ' 12"E  Longitude 


Output 

Data/Units 


2  MOO 

49400 

42.44S7 

-41.07)2 

-27.0007 

102.2712 
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Program  Listings 


Master  Data  Card  Program 


000 

76 

LBL 

050 

75 

100 

55 

Oul 

42 

7  TO 

051 

02 

101 

02 

2 

0  0 2 

8T  * 

052 

0  9 

9 

102 

95 

= 

0  0 

00 

00 

053 

os 

r. 

103 

42 

STD 

004 

97 

DS2 

054 

Q 

104 

01 

01 

-  “  LT 

f  J  i  J 

0  0 

00 

055 

02 

1  05 

4  3 

RCL 

0  0  0 

42 

a 
1 — 

056 

06 

6 

106 

24 

24 

007 

0  0 

o 

057 

*i  er 

l/X 

107 

75 

- 

ooo 

4!  4! 
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Master  Data  Card  Program  (continued) 
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Operational  Data  Card  Program 
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Operational  Data  Program  Card  (continued) 
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Position  Determination  Program  (continued) 
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ion  Program  (continued) 
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D.  Loran-C  Fixing  Algorithms 

The  development  of  the  Loran  fixing  algorithms  in  this 
report  is  presented  in  more  detail  in  a  companion  report 
{Ref.  3]  and  will  not  be  repeated  here. 

The  basic  Loran-C  equation  [Ref.  4]  can  be  written  as 

T  =  [Tg  +  p(Tg)l  -  [Tm  +  p(TMn  +  [Tb  +  p(TN)]  +  6  (1) 

where 

T  is  the  "indicated  time  difference”  in  microseconds, 
T„,T,,  is  the  distance,  in  microseconds,  from  the  master 

M.  b 

and  the  slave  to  the  receiver,  respectively, 

Tb  is  the  distance,  in  microseconds,  between  the 
master  and  the  slave, 

<5  is  the  assigned  coding  delay,  in  microseconds,  and 
p (T)  is  the  secondary  phase  correction,  in  microseconds, 
for  an  all  sea  water  path  of  length  T. 

The  quantity 

At  =  [Tg  +  p(TB)]  +  6 

is  a  constant  for  each  master/slave  pair.  The  following 
World  Geodetic  System  1972  (WGS  72)  values  have  been  adopted 
for  Loran-C  navigation  {Ref.  4]: 
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vQ  =  299792458  meters/second  is  the  velocity  of  light 
in  free  space, 

ii  =  1.000)38  is  the  index  of  refraction  of  the  surface 
of  the  earth  for  standard  atmosphere  and  100kHz 
elect  'omayentic  waves, 

a^  =  6378135.00  meters  is  the  equatorial  radius  of  the 
earth,  and 

f  =  1/298.26  is  the  flattening  factor  (1  -  b/ae,  where 
b  is  the  polar  radius)  of  the  earth. 

Accurate  formulas  for  computing  the  secondary  phase 
correction  p(T)  are  contained  in  Reference  4,  but  for  use 
with  the  handheld  calculator  the  following  linear  approximation 
[Ref.  3]  will  be  used: 


where 

and 


P (T)  =  a1  +  a2T  , 


a^  =  -0.321, 
a2  =  0.000635. 


Using  this  approximation,  it  is  possible  to  solve  Equation  1 
for  the  quantity  Tg  -  T^.  We  find 


Ts  -  TM  =  (T  -  At ) / ( 1  +  a2) 


(2) 


On  the  surface  of  a  sphere  a  hyperbolic  line  of  position 


can  be  represented  by  the  equation  [Ref.  3,  page  175] 


tan  r 


(3) 


_ cos  2a  -  cos  2c 

sin  2c  cos  u  +  c  sin  2a 

where  the  origin  of  the  coordinate  system  is  at  the  prime  focus 
of  the  spherical  hyperbola,  2c  is  the  spherical  arc  joining 
the  foci,  2a  is  a  constant  for  any  one  LOP,  and  r  and  w 
are  the  spherical  coordinates  of  a  point  on  the  LOP.  If  the 
base  line  of  the  coordinate  system  is  the  arc  joining  the  foci 
then  co  is  the  spherical  polar  angle  from  the  base  line  to  a 
point  P  on  the  LOP  and  r  is  the  spherical  polar  distance 
(or  arc)  from  the  prime  focus  to  P.  Using  the  Loran  system 
we  take  c,  =  +1  if  the  prime  focus  is  at  a  master  station 
and  we  take  r,  =  -1  if  the  prime  focus  is  at  a  slave  station. 

If  we  take  v  =  vQ/n  to  be  the  velocity  of  100kHz  electro 
magnetic  radiation  of  the  earth's  surface  then 

2a  .  v (Tg  -  TM)/ae  , 

or,  using  Eq .  (2) , 

2a  -  (T  -  At) /a  ,  (4) 

where 

a  (1  +  a2) 

a  =  - 7 — ±-  =  21295.87  ps  . 

P  vQ/n 

The  baseline  between  master  and  slave  can  be  obtained  from 

2c  =  v  T_/a  .  (5) 

O  6 

Here  2c  is  computed  by  program  card  1  (preparation  of  master 
data  cards)  using  the  algorithm  in  Section  E. 
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Consider  a  Loran-C  triplet  with  master  stations  colocated. 
Let  •'  ^  and  2  denote  the  azimuth  angles  of  slave  1  (S^) 
and  slave  2  (S^) ,  respectively,  measured  from  North  toward  the 
East  from  the  master  stations  (M)  (see  Fig.  3).  Further,  let 
it  and  r  denote  the  azimuth  and  spherical  polar  arc  (distance) 
of  the  receiver  (R)  from  M.  For  this  geometry,  Eq.  (3)  can 
be  written  as 


where 


tan  r^ 


a . 

_  x  _ 

C~  cosTa-T~l  +—A~ 


(6) 


and 


> 

II 

c, . 

X 

sin 

Bi  - 

cos 

2a . 

l 

II 

•H 

u 

sin 

2c  . 
l 

.  th 


for  the  i —  Loran  pair,  i  =  1,2 
can  be  eliminated  in  Eq.  (6). 
rewritten  as 


.  Since  r  =  r^  *  r2,  tan 
The  resulting  equation  can  be 


C  cos  u  +  S  sin  a  =  K  , 


(7) 


where 


C 

-  Blc2 

cos  ^2 

-  b2Ci 

cos  t, 

S 

=  3^2 

sin  1 2 

B2C1 

sin  C 

and 

K 

= 

-  B1A2 

• 
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If  we  define  >  0  and  t  by  the  equations 


i  cos  i  -  C, 


then 


i  sin  i  =  S, 


n  =  /c2  +  S2  , 


r  =  qatn(S,C) 


Here  the  function  gatn(y,x)  is  the  arctangent  of  y/x 
adjusted  for  the  proper  quadrant  according  to  the  signs  of  x 
and  y.  /A  compact  form  of  this  function  is 


where 


qatn(yfx)  =  tan - * -  +  nt(x  <  0?) 

x  +  10~  t  (x  -  0?) 


t(z)  =  1  when  z  is  true 


t  ( z)  =  0  when  z  is  false. 


When  convenient  we  will  use  the  notation  qatn(y/x)  inter¬ 
changeably  with  qatn(y,x) .  The  qatn  function  is  equivalent 
to  the  polar  angle  obtained  using  the  rectangular  to  polar 
conversion  function  on  the  TI-59. 

Now  substitute  Eq.  (8)  into  Eq.  (7)  and  solve  for 


fc; 


a  =  )  +  cos  1 ( k /p ) 


1 


to  obtain  the  azimuth  angle  <*  of  the  two  points  of  intersection 
of  the  LOP's.  finally  we  obtain  a  value  for  r  by  substituting 
each  “  into  Eg.  (5).  We  find  that 


r 


qatn 


B. 

l 


cos 


for  i 


1  or  2 . 


The  distance  and  azimuth  from  M  or  the  triplet  vertex  can  be 
converted  into  the  latitude  and  longitude  of  the  two  possible 
positions  of  R. 

The  fixing  algorithm  then  uses  a  and  r  in  the  direct 
solution  algorithm  of  spheroidal  geodesy  (Section  F) . 
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E.  The  Rt'Vi  i  Sf  (  I  nvi  rst.' )  Solution  A  1  gor  l  thin 

This  ’■  solution  algorithm  is  u  modi  t  nratiun  of  the 

first  order  in  tl.t  toning  (t)  algorithm  yiven  by  Thomas  [Red.  5, 
pp.  8-10).  Thomas'  notation  has  been  followed  as  closely  as 
possible  for  ease  of  comparison  of  the  algorithms.  The  yatn 
function  is  defined  in  Section  D.  West  longitudes  (')  and 
South  latitudes  (|)  are  negative.  We  are  given  the  points 
P^(4*^#A^),  ?2^2,A2^  °n  t^ie  slJf*er°id  and  are  to  find  the 
distance  S  between  the  points  and  the  forward  and  back 
azimuths,  and  ll21’  Given  quantities  are  <|ij,  A^, 

and  A 2 •  No  assumptions  about  the  relative  location  of  P^ 
and  P2  are  required.  The  modified  I'everne  solution  algorithm 
is : 

tb  =  tan  *1(1-0  tai  4ml,  i  =  1,  2, 

=  («.  +  0  )/2,  At  =  (0  -0,)/2,  AA  =  A  -  A., 
m  1  2  m21  21 

2  2  2  2 

AA  *  AA/2,  H  =  cos  A0  -  sin  8  =  cos  0  -  sin  A0  =  cos  0,  cos  8  , 

m  m  m  m  m  1  2 

L  *  sin2  A0  +  H  sin2  AA  =  sin2(d/2),  1  -  L  ■  cos2(d/2), 
m  m 

-1  2  2 
d  *  cos  (1  -  2L) ,  U  *  2  sin  0  cos  A0  /(I  -  L) , 

m  01 

.  2  2 

V  »  2  sin  A0  cos  0  /L,  X  *  U  +  V,  Y  =  U  -  V, 
w  m 

T  -  d/sin  d,  6 jd  -  f (TX  -  Y)/4,  S  *  ag (T-6 jd)  sin  d, 

F  -  2(Y  -  (1  -  21.)  (4  -X)],  G  *  f T/2 , 

Q  =  -(FG  tan  Aa)/4,  AA '  =  (AA  +  Q) / 2 
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t.  =  uutn (-sin  AO  cos  AX',  cos  6  sin  AX'), 

1  m  m  m  m 

t  =  out n (cos  AO  cos  Aa  1  ,  sin  o  sin  AX')  , 

2  1  rn  in  m  m 

“l2  =  fcl  +  V  “21  =  S  ’  C2  * 

This  reverse  solution  algorithm  is  used  by  program  card  1 
(preparation  of  master  data  cards)  to  compute  the  baseline 
distance  2c  and  the  azimuths  f^g  and  between  the 

master  and  slave  stations  of  a  Loran  pair. 

Details  of  the  modifications  made  to  Thomas’  algorithm 
are  contained  in  Reference  3. 


i 


i 


( 


i 


* 
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F.  The  Direct  Solution  Algorithm 

This  solution  algorithm  is  a  modification  of 

the  first  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  7-8).  Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 

The  gatn  function  is  defined  in  Section  D.  West  longitudes 
and  South  latitudes  are  negative.  We  are  given  the  point 
Pl<^l'^l)  on  the  sPherio:i-d*  where  4^,  *  are  the  geodetic 
latitude  and  longtiude  (geographic  coordinates) ;  the  forward 
azimuth  u12  and  the  distance  S  to  a  second  point  p 2  < ^' 2  '  ' 2 ^  '* 

and  from  these  we  are  to  find  the  geographic  coordinates 

^2'  ^2  anC*  t*ie  azimuth  '*21  •  T^e  9iven  quantities  are 

^1'  *i#  ai2  and  S‘  N°  assurnPt^ons  about  the  relative  loca¬ 
tion  of  P1  and  P2  are  required.  The  modified  ,!ira-i 

solution  algorithm  is: 

=  tan  1 1  ( 1  - f )  tan  4^1.  M  =  cos  ti^  sin  uJ2 

N  *  cos  cos  u12'  ci  =  fM'  c2  =  ~  m2,/,4> 

D  =  1  -  2c 2  -  c ^ M#  P  =  c2/D,  o  =  qatn(N,  sin  0^ 

d  *  S/(a  D)  ,  u  *  2(0, -d)  ,  W  =  1  -  2P  cos  u, 
e  1 

V  =  cos(u  ♦  d) ,  Y  =  2PVW  sin  d,  A a  =  d-Y, 

{X2\  “  qatnl-M,  -  (N  cos  Ao  -  sin  sin  Ao ) 1  , 

2  2  1/2 
K  =  (1-f) (M  +  (N  cos  Ao  -  sin  sin  Ao )  ]  , 

4  ?  =  tan  ^IlsinO^  cos  A o  +  N  s i n  Ao  )  /K 1  , 

Ari  =  qatn<sin  Ao  sin  cos  0  {  cos  Ao  -  sin  0^  sin  Ao  t.nti  u  ), 

«i  -  CjAo,  AA  *  An  -  H,  A2  =  +  AA. 
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This  direct  solution  algorithm  is  used  by  program 
card  3  (improved  fix  program)  to  compute  the  latitude  and 


longitude  of  the  receiver  using  the  azimuth  and  range  of  the 
receiver  from  the  Loran  triplet  vertex. 

Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3. 


G.  Some  Results 


The  following  procedure  was  used  to  test  the  fixing 
algorithm:  A  more  accurate  version  of  the  reverse  solution 

algorithm  (see  Reference  3)  was  used  to  compute  time  delays 
at  selected  positions  from  Loran-C  pairs.  These  computed 
time  delays  do  not  depend  upon  a  fixing  algorithm.  The 
fixing  algorithm  was  tested  by  determining  the  distance 
between  the  fix  and  the  positions  used  to  compute  the  time 
delays.  Typical  results  are  presented  in  Tables  I  through 
IV.  As  can  be  seen,  all  differences  (A)  are  within  1  n.mi. 
The  greatest  difference  (A  =  0.92  n.mi.  in  Table  IV)  was 
obtained  at  a  pdsition  located  2270  n.mi.  from  the  triplet 
vertex. 


Table  I.  9940X/9940Y 


Position 

Lat  Long 

Calculated 

9940X 

Time  Delay 

9940Y 

TI-59 

Lat  (N) 

Fix 

Long  (W) 

A 

n.mi 

24°N 

122°W 

27726.19 

40912.76 

23°59 ' 55" 

122°00 ' 01" 

0.08 

26 

122 

27715.97 

40998.39 

25°59,57" 

122°00'01" 

0.05 

28 

122 

27702.41 

41117.84 

28°00 1 00" 

122°00 1 00" 

0.00 

30 

122 

27683.53 

41291.85 

29°59 1 59" 

122°00 1 00" 

0.02 

32 

122 

27655.47 

41555.46 

32°00 1 00" 

122°00'00" 

0.00 

34 

122 

27609.63 

41959.57 

34"00'00" 

122°00 ’00" 

0.00 

36 

122 

27523.56 

42544.11 

seooo'oo" 

121°59 ' 59" 

0.01 

38 

122 

27334.61 

43248.22 

38°00 ' 00" 

122°59'59" 

0.01 

Table  II.  9940Y/9940W 


Position 

Calculated 

Time  Delay 

TI-59 

Fix 

A 

n.mi 

Lat 

Long 

9940Y 

9940W 

Lat  (N) 

Long  (W) 

37°N 

122°W 

42892.86 

16257.23 

36°59 ' 59" 

122°01 ’01" 

0.02 

37 

125 

43056.68 

15765.13 

37°00 ' 00" 

125"00'00" 

0.00 

37 

128 

43137.78 

15327.12 

37°00'00" 

128°00 ' 00" 

0.00 

37 

131 

43191.10 

14970.77 

37°00 ' 00" 

131°00'00" 

0.00 

37 

134 

43232.38 

14683.74 

O 

O 

O 

o 

o 

134°00 '00" 

0.00 

37 

137 

43267.42 

14449.40 

37°00 ' 00" 

137°00 1 00" 

0.00 

37 

140 

43298.80 

14254.02 

37°00,00" 

140°00 ' 00" 

0.01 

37 

143 

43327.85 

14087.43 

37°00‘01" 

142°59 1 59" 

0.02 

Table  III.  7930Z/9930X* 


Position 

bat  Long 

Calculated 

79302) 

Time  Delay 

9930X 

TI-59 

Lat  (N) 

Fix 

Long  (W) 

n  .mi 

60°N 

30°W 

52437.86 

28451.72 

60°00  1 03" 

29°59' 32" 

0.24 

58 

35 

51960.93 

28391.50 

58°00 1 00" 

34  ’’59  1 48" 

0.09 

56 

40 

50992.37 

28359.15 

39°59 1 54" 

0.06 

54 

45 

49292.46 

28370.85 

53°59'59" 

44° 59 1 57" 

0.03 

52 

50 

47165.60 

28490.64 

52  °00 ' 00" 

49° 59 • 57" 

0.01 

50 

55 

45236.59 

29070.48 

50°00 ' 00" 

55  °00 ' 00" 

48 

60 

44505.60 

30991.94 

O 

o 

o 

o 

0 

CO 

60  °00 ' 00" 

46 

65 

44475.70 

33697.14 

46°00 1 00" 

03 

U1 

o 

O 

o 

o 

o 

44 

70 

44588.91 

36567.42 

43°59 1 59" 

69  °59 ' 59" 

Table  IV.  9930W/9930X* 


Position 

Lat  Long 

Calculated 

9930W 

Time  Delay 

9930X 

TI-59 

Lat  (N) 

Fix 

Long  (W) 

A 

n  .mi 

9®N 

47°W 

13058.04 

36466.46 

8°59 ' 59" 

46°59 

22" 

0.92 

12 

52 

12984.71 

37288.35 

11 °59 ' 34" 

51°59 

37" 

0.57 

15 

57 

12898.73 

38267.58 

14°59'45" 

56°59 1 47" 

0.28 

18 

62 

12793.91 

39431.32 

17°59 ' 52" 

61  °59 

54" 

0.16 

21 

67 

12656.52 

40794.36 

20°59 ' 56" 

66°59 

57" 

0.08 

24 

72 

12451.30 

42330.55 

23°59 ' 59" 

71  ”59 

59" 

0.02 

27 

77 

12097.12 

43876.62 

27°00 ' 01" 

77  °00 

00" 

0.02 

30 

82 

12973.95 

44768.53 

30°00’01" 

82  °00 

06" 

0.09 

*The  9930  system  was  disestablished  30  September  1979.  It  is  included  here 
for  illustration  only. 


36 


H.  References 


1.  J.  A.  Pierce,  A.  A.  McKenzie,  and  R.  H.  Woodward,  editors, 
loran ,  M.I.T.  Radiation  Laboratory  Series,  McGraw-Hill 
Book  Company,  Inc.,  1948. 

2.  G.  Hefley,  The  Development  of  Loran-C  Navigation  and  Timing, 
National  Bureau  of  Standards  Monograph  129,  U.  S.  Depart¬ 
ment  of  Commerce,  U.  S.  Government  Printing  Office, 
Washington,  D.  C.  20402,  October  1972. 

3.  R.  H.  Shudde,  "An  Algorithm  for  Position  Determination 
Using  Loran-C  Triplets  with  a  BASIC  Program  for  the 
Commodore  2001  Microcomputer, "  Technical  Report 
NPS55-80-009 ,  March  1980,  Naval  Postgraduate  School, 
Monterey,  CA  93940. 

4.  LORAN  HYPERBOLIC  LOP  FORMULAS  AND  GENERAL  SPECIFICATIONS 
for  loran-c  (20  June  1949)  were  obtained  from  G.  R.  Young, 
Acting  Chief,  Navigation  Department,  Defense  Mapping 
Agency,  Hydrographic/Topographic  Center,  Washington,  D.C. 
by  private  communication,  5  March  1980. 

5.  Paul  D.  Thomas,  "Spheroidal  Geodesics,  Reference  Systems, 
and  Local  Geometry,"  PS-138,  U.  S.  Naval  Oceanographic 
Office,  Washington,  D.C.,  January  1970. 

6.  R.  H.  Shudde,  "Position  Determination  with  Loran-C 
Triplets  and  the  Hewlett-Packard  HP-67/97  Programmable 
Calculators,"  Technical  Report  NPS55-80-010 ,  March  1980, 
Naval  Postgraduate  School,  Monterey,  CA  93940. 


37 


APPENDIX.  Loran-C  Station  Parameters 


The  following  list  contains  the  pertinent  parameters 
for  each  Loran-C  station  pair.  This  list  was  compiled  from 
data  in  Reference  4.  Each  row  contains  the  following  inform¬ 
ation: 

1.  The  Loran-C  station  pair  designator 

2.  At,  the  sum  of  the  coding  delay  plus  one  way  base¬ 
line  time,  including  the  secondary  phase  correction 
for  an  all  seawater  path,  in  micrseconds. 

3.  The  master  station  latitude. 

4*  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  report,  the  convention  that  a  minus  (-)  sign  denotes 
South  latitudes  and  West  longitudes  has  been  followed 
throughout  Some  users  may  wish  to  reverse  the  longitude 
convention  so  that  (-)  longitudes  are  East  rather  than  West. 
This  may  be  accomplished  by  consistently  changing  the  sign 
of  the  longitude  everywhere  that  it  is  required  for  data 
input;  no  program  modifications  are  required.  If  the  sign 
convention  for  longitude  is  reversed,  care  should  be  taken 
that  the  signs  of  all  longitudes  in  the  list  are  reversed. 

In  columns  3  through  6  the  latitudes  and  longitudes  appear 
to  be  in  decimal  form,  but  the  actual  format  is  DDD.MMSSFF 
where 
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DDD  designates  degrees, 

MM  designates  minutes, 

SS  designates  seconds,  and 

FF  designates  hundredths  of  seconds 
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Coverage  of  Loran-C  Systems 


Station 

Location 

4990 

Central  Pacific 

59  30 

East  Coast,  Canada 

5990 

West  Coast,  Canada 

7930 

North  Atlantic 

7960 

Gulf  of  Alaska 

7970 

Norwegian  Sea 

7980 

Southeast  U.S.A. 

7990 

Mediterranean  Sea 

8970 

Great  Lakes 

9930 

East  Coast,  U.S.A. 

9940 

West  Coast,  U.S.A. 

9960 

Northeast  U.S.A. 

9970 

Northwest  Pacific 

9990 

North  Pacific 
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